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We discuss the low-temperature dynamics of magnetic flux lines in samples with a family of parallel twin planes. A current 
applied along the twin planes drives flux motion in the direction transverse to the planes and acts like an electric field applied 
to one-dimensional carriers in disordered semiconductors. As in flux arrays with columnar pins, there is a regime where the 
dynamics is dominated by superkink excitations that correspond to Mott variable range hopping (VRH) of carriers. In one 
dimension, however, rare events, such as large regions void of twin planes, can impede VRH and dominate transport in samples 
that are sufficiently long in the direction of flux motion. In short samples rare regions can be responsible for mesoscopic effects. 
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The static and dynamic properties of magnetic flux 
lines in copper-oxide superconductors are strongly af- 
fected by pinning by point, linear and planar disorder 
0. Linear correlated disorder in the form of columnar 
defects produced by bombardment of the crystal with 
energetic heavy ions has been shown to greatly enhance 
pinning in both yttrium and thallium-based compounds 
0. Twin boundaries are an example of planar disor- 
der that is ubiquitous in superconducting Y Bo2CuzOt- x 
and LaiCuO/±. Early decoration experiments indicated 
that the superconducting order parameter is suppressed 
at a twin boundary and the twin attracts the vortices || . 
Extensive investigations of twin-boundary pinning have 
been carried out by Kwok and coworkers ||||. These 
authors studied a variety of YBCO single crystal sam- 
ples containing single families of parallel twins lying in 
planes spanned by the c axis, with spacings ranging from 
microns down to several hundred Angstroms. Trans- 
port experiments [Q show clear evidence of strong twin- 
boundary pinning even in the flux liquid phase for ex- 
ternal fields along the c axis and driving currents in the 
ab plane and parallel to the plane of the twins (resulting 
in a Lorentz force normal to the twin planes). For this 
geometry the linear resistivity drops sharply at a char- 
acteristic temperature where twin-boundary pinning sets 
in Q]. In addition, there is a sharp downward dip in the 
resistivity as a function of angle as the external field is 
rotated through the c direction [||. 

The static and dynamic properties of flux-line assem- 
blies in the presence of a random array of columnar pins 
have been studied by mapping the physics of magnetic 
flux lines onto the problem of localization of quantum 
mechanical bosons in two dimensions At low tem- 
peratures there is a "Bose glass" phase, with flux lines 
localized on columnar pins, separated by a phase transi- 
tion from an entangled flux liquid of delocalized lines. 



Transport in the Bose glass phase closely resembles the 
variable-range hopping (VRH) of electrons in disordered 
semiconductors in two dimensions 

In this paper we study flux-line dynamics in the pres- 
ence of a single family of parallel twin boundaries lying in 
planes containing the c axis, for H || c. We find that due 
to the quasi-one-dimensional nature of vortex transport 
a new regime can arise at low currents where flux-line 
dynamics is dominated by rare events, such as large re- 
gions voids of twin planes, that can be responsible for 
new mesoscopic phenomena. 

At low temperatures, when the average vortex spac- 
ing ao « (<po/B) exceeds the average distance d between 
twin planes, all flux lines are localized by the pinning po- 
tential in the direction normal to the twins, progressively 
"filling" the planar pins as the field is increased. We only 
consider fields below Bf w (f>o/d 2 , where the flux lines fill 
the twins, and neglect any additional weak point disorder 
in the sample. We focus on flux motion transverse to the 
twin planes, which resembles the hopping of electrons in 
one- dimensional disordered superconductors. The cur- 
rent density in the usual hopping conductivity problem 
corresponds to the vortex velocity (i.e., voltage) and the 
electrical conductivity maps onto the resistivity from vor- 
tex motion. Transport in this low temperature and field 
regime is dominated by single- vortex dynamics. Flux mo- 
tion is described in terms of thermally activated jumps 
of the vortices over the relevant pinning energy barriers 
U(L, J), yielding a resistivity p = 8/J, 

p(T) « p,e- u ' T , (I) 

where po is a characteristic flux-flow resistivity. In the 
following we determine the barrier heights U(L,J) (Table 
I) corresponding to various transport regimes and the 
boundaries between the different regimes in the (L, J) 
plane. 
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Our starting point is the model-free energy discussed 
in H for flux lines defined by their trajectories {r^z)} 
as they traverse a sample of thickness L with a family of 
parallel twin boundaries parallel to the zx plane, 
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where V(r) is the intervortex interaction 
single-flux line free energy, 
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Here ei is the tilt modulus (?i = (M±/M z )ei < ei, 
with ei the line tension) and Vd (y) represents a z-and x- 
independent pinning potential At low temperatures 
we model Vd (y) as an array of identical one-dimensional 
square potential wells of depth l/ , width 2b and average 
spacing d >> bo, passing completely through the sample 
in the x and z directions and centered at uniformly dis- 
tributed random positions along the y axis. 

As discussed in Q and JhJ, many relevant results re- 
garding the statistical mechanics of flux lines can be ob- 
tained from elementary quantum mechanics by mapping 
the vortex trajectories onto the imaginary time path in- 
tegrals of two-dimensional particles in a static random 
potential Vd(v)- In this mapping fc^T plays the role of h, 
ii that of the mass of the fictitious quantum particle, L _1 
that of the particle's temperature. If we assume that each 
flux line spends most of its time near the attractive twin 
planes, the dynamics of flux lines driven by a Lorentz 
force normal to the twin planes can be described by a 
tight-binding model for one- dimensional bosons, where 
the lattice sites are the positions of the twin planes. Flux 
lines can "tunnel" between twins i and j separated by a 
distance dij at a rate given by the tunneling matrix ele- 
ment Uj ~ 2U(T)e- E *^ T , where E tj = ■ y /2e 1 U(T)d ij is 
the energy of a "kink" configuration connecting the two 
pins and U(T) is the effective binding free energy per unit 
length of a flux line trapped near a twin plane. At low 
temperature U(T) « Uq. Thermal fluctuations renormal- 
ize Uq above a crossover temperature, as discussed in [ffo|. 
The modeling of the dynamics of two-dimensional bosons 
in terms of a tight binding model in one dimension can 
be understood as a result of vortex-vortex interactions 
which confine each flux line within a "cage" of radius 
~ a formed by its neighbors. Each line can be viewed 

TABLE I. Energy barriers determining the various contri- 
butions to the resistivity of Eq. 1, with a = g(jj,)dU(T). 



Linear 


Nonlinear 


U rf = UL = E k (L/w k ) 


Uhi = S fc (Ji/J) 


U-nnh = Ek 




U Mo tt = E k (L/aw k ) 1/2 


Uvrh = E k (Ji/aJ) 1/2 



as moving within a one-dimensional "channel" of width 
~ oq. Another effect of interactions is the energy cost for 
an additional flux line to be placed on an already filled 
twin boundary. This is incorporated into an energy cost 
Vb for double occupancy of a site of the one-dimensional 
tight-binding lattice, which is estimated as Vb ~ 4eod 2 ja\ 
for T < T*d/b [ln(a /2d)]- 1 and d « a « X ab . 

We consider vortex transport in the presence of a driv- 
ing current J _L H parallel to the twin planes, i.e., 
J = — Jx. The applied current exerts a Lorentz force 
per unit length on the vortices, = — z x J = yf L in 
the direction transverse to the twin planes. In the context 
of the analogy with boson quantum mechanics, this term 
represents a fictitious "electric field" E = ^z x J = yJ/c 
acting on particles with "charge" <fio- At low temper- 
atures the critical current can be obtained by equat- 
ing the Lorentz force to Uo/bo, where bo ~ £ a 6, i.e., 
J c (0) w cUo/cf>obo. Renormalization of J c by thermal 
fluctuations have been discussed in |l. 
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FIG. 1. The (L, J) phase diagram for a — g(fx)dU < 1. 

If we neglect rare events, the energy barriers deter- 
mining the resistivity of Eq. [I] can be estimated as the 
saddle-point free energy associated with the low-lying ex- 
citations from the ground state (where all the lines are 
localized on twins). The results are summarized in Ta- 
ble 1. The phase diagram in the (L, J) plane is con- 
trolled by the parameter a = g((i)dU(T), where g(n) 
is the density of states for the most weakly bound flux 
lines and fj, ~ <f>o{H — H c i)/Ait is the chemical potential 
which fixes the flux line density. The low-lying excita- 
tions that govern transport at low temperature have been 
discussed before in the context of flux arrays pinned by 
columnar defects ||. The only difference here is in the 
exponent of the barrier associated with VRH, both in 
the linear and nonlinear regime. In the presence of lin- 
ear defects Uvrh ~ (J/J) 1 / 3 , while for planar defect 
the corresponding exponent is 1/2, as shown in Table 1. 
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The boundaries of the (L, J) phase diagrams (Figs. 1 
and 2) have not been discussed before in detail. The line 
J = Jl(L) defines the boundary in the (L, J) plane that 
separates the regions of linear (J < Jl(L)) and nonlin- 
ear (J > Jl(L)) response, with Jl — cEk / '(foLd) and 
Ek = \fTjJ(T)d the energy of a kink connecting neigh- 
boring pins separated by d. In the thermodynamic limit 
Jl — * and the IV characteristic is nonlinear at all cur- 
rents. In samples of very small thickness L there is a 
linear resistivity at low currents due to the flow of flux- 
line segments of length L and typical transverse width 
y r f « dL/wk, with Wk = E^/U the width of a kink 
connecting pins separated by d. When < L < L±, 
where L\ = E^/j is the length below which dispersion 
from tunneling and interactions can be neglected, trans- 
port occurs via the hopping of vortices between near- 
est neighbor (nn) pinning sites. Here 7 is the width 
of the impurity band arising from level dispersion from 
tunneling and intervortex interaction. The region of the 
phase diagram dominated by nn hopping is only present 
if L\ > Wk, or 7 < U (ll). When L > L\ the disper- 
sion of energies between different pinning sites makes mo- 
tion by nearest neighbor hopping energetically unfavor- 
able. Tunneling occurs instead via the formation of "su- 
perkinks" that throw a vortex segment onto a spatially 
remote pin connecting states which optimize the tunnel- 
ing probability. Tunneling via supcrkinks is the analogue 
of Mott variable-range hopping of electrons between lo- 
calized states in semiconductors. The corresponding en- 
ergy barriers Umou is given in Table 1. 

For J > Jl(L) the resistivity is nonlinear. For J\ < 
J < J c , with Ji = cU(T)/(f)od, flux motion occurs via 
thermally activated "half-loop" configurations of trans- 
verse width yhi ~ U/ II- For J < J\ the size of the 
transverse displacement of the liberated vortex segment 
exceeds the average distance d between twin planes and 
transport occurs via VRH which generalizes the Mott 
mechanism to the nonlinear case. The characteristic cur- 
rent scales that governs VRH is Jo = J\/a. The VRH 
contribution to the resistivity dominates that from half 
loop only if if Uvrh < Um, or J < J2 = J\u. All current 
scales in Fig. 1 are much smaller than the pair break- 
ing current J p b = 4ceo/(3v / 30o£a&)- The Mott and the 
rigid flow regimes are separated by a horizontal line above 
which Umou < U r f. Similarly, the condition Uvrh = Uhi 
yields the vertical line separating the VRH and half-loop 
regions. 

Transport in flux arrays with columnar pins is de- 
scribed by a phase diagrams similar to that of Fig. 1. 
The difference for planar disorder is that transport is 
one dimensional in this case and rare fluctuations in the 
spatial distribution of twins can impede VRH and domi- 
nate flux-line dynamics. The vortex line can encounter a 
region where no favorable twins are available at the dis- 
tance of the optimal jump. The vortex will then remain 
trapped in this region for a long time and the resistivity 



can be greatly suppressed. Rare fluctuations can also oc- 
cur in samples with columnar pins, but in that case they 
will dominate transport only at extremely small fields, 
when the number of rare regions exceeds the number of 
vortices. 

At a given temperature and for applied currents below 
Jl, a vortex can jump from one twin plane to another 
at a distance y if the energy difference per unit length 
between the initial and final configuration is within a 
range Ae ~ Eky/Ld. A trap is then a region of con- 
figuration space (y, e) void of localized states within a 
spatial distance y and an energy band Ae around the 
initial vortex state. A vortex that has entered such 
a trap or "break" will remain in the trap for a time 
t w w to exp(2y/7j_), where l± = T/\f2c\U is the localiza- 
tion length and to a microscopic time scale. The proba- 
bility of finding such a break is given by a Poisson dis- 
tribution, P(y) ps P (y)exp[-Ag(iJ 1 )yAe), where P (y) is 
the concentration of localized states in the energy band 
Ae, Po(y) ~ 2Ag(/j,)Ae and A ~ 1 a numerical constant. 
The mean waiting time between jumps is 

/>oo 

t^w / dyP(y)t e 2 ^ T \ (4) 
Jo 

For L » (T/Ek) 2 awk = L* , the integral can be evalu- 
ated at the saddle point, corresponding to the situation 
where the mean waiting time is controlled by "optimal 
breaks" of transverse witdth y\ sa l±L/L*, with the re- 
sult, t w ~ fo\/ L/ L*e L / L . The optimal breaks corre- 
spond to the longest trapping time and are most effective 
at preventing flux motion. The inverse of the trapping 
time determines the characteristic rate of jumps, yielding 
a linear resistivity, 

p W ^Po(T/W Uoi t)e- ( """" /T)2 . (5) 

For currents above Jl the typical energy per unit 
length available to a flux line for jumping a distance y 
is Ae ~ f L y. Again, for J « (J^a^Ek/T) 2 = J* 
the main contribution to the resistivity is from "optimal 
traps" of transverse size yl ~ l±(J* /J), with the result, 

p b pb po(T/U V RH)e- {UvRH/T)2 . (6) 

The contribution to the resistivity from tunneling a la 
Mott (both in the linear and nonlinear regimes) always 
dominates that from hopping between rare optimal traps 
if both mechanisms of transport can occur. On the other 
hand, in one dimension if the sample is wide enough in 
the direction of flux-line motion to contain optimal traps, 
tunneling a la Mott simply cannot take place because flux 
lines cannot get around the traps. These rare traps with 
large waiting times will then control transport. If W is 
the sample width in the y direction, the condition for 
having optimal traps of width y\ h is P(y1 b )W > 1. 
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FIG. 2. The (L, J) phase diagram for a = g(p)dU < E k /T. 

Optimal traps will be present only if J > J w — 
J* I ln(2W/l ± ) for J> J L and if L < L w = L* ln(2W/7j_) 
for J < J L . 

A flux line can, however, escape a break by nucleating 
a half-loop, i.e., by tunneling directly into conduction 
band, if yu > y%, or a < Ek/T. The only effect of 
rare fluctuations in this case is of suppressing VRH for 
J > J w , extending to lower currents the region where 
transport occurs via half-loop nucleation. For instance 
the right (high current) boundary of the VRH region 
in the phase diagram of Fig. 1 will be pushed down 
to J w if J w < J 2 , or ln(2W/y T ) > (E k /aT) 2 . Simi- 
lar considerations apply to the linear response. Only if 
a > Ek/T, there will be a portion of the (J, L) plane 
where breaks dominate transport, as shown in Fig. 2. 
For YBCO, we estimate Ek ~ lKA~ 1 d. Assuming 
a ~ C//7, the condition a > Ek/T can only be satis- 
fied at low fields (B < IKG for d ~ 200l). The sample 
will contain optimal breaks if W > 30Aexp(J*/J), with 
J*~4x 10 5 Amp/cm 2 at 80K. 

If the sample is too short to contain optimal breaks, 
i.e., WP(y ~ y'i b ) < 1, the dynamics is controlled by the 
trap with the longest waiting time, ~ exp(j/f /l±), 

with yj determined by the condition WP(yf) ~ 1. The 
corresponding resistivity is pw ~ Po e ~ y f' lx . In this case 
the relevant physical quantity is the logarithm of the re- 
sistivity, 

Mpw/po) = ~Vf/l± 

The leading dependence of Eq. |?j on current and tempera- 
ture is the same as that of the VRH contribution. Equa- 
tion ([?]) also contains, however, logarithmic terms that 
in sufficiently short samples will give a random spread of 
values of the resistivity from sample to sample. These 
effects have been discussed for semiconductors fl3l. In 



this case a more relevant physical quantity rather than 
the resistivity itself is the distribution of the logarithms 
of the resistivity over different samples. The expression 
determines the position of the maximum of this dis- 
tribution. 
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